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Abstract 



A search is presented for free heavy long-lived fractionally charged particles pro- 
duced in pp collisions at \/s = 7TeV. The data sample was recorded by the CMS 
detector at the LHC and corresponds to an integrated luminosity of 5.0 fb -1 . Candi- 
date fractionally charged particles are identified by selecting tracks with associated 
low charge measurements in the silicon tracking detector. Observations are found to 
be consistent with expectations for background processes. The results of the search 
are used to set upper limits on the cross section for pair production of fractionally 
charged, massive spin-1/2 particles that are neutral under SU(3)c and SU(2)l- We 
exclude at 95% confidence level such particles with electric charge ±2e/3 with masses 
below 390 GeV, and those with charge ±e/3 with masses below 250 GeV. 
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The reasons for expecting physics beyond the standard model to manifest itself in pp collisions 
at the Large Hadron Collider (LHC) are as compelling as ever. As suggested in [1], one example 
of physics beyond the standard model that may have eluded previous searches is a new particle 
with an electric charge less than that of the electron. Owing to their lower ionization energy 
loss, the trajectories of such fractionally charged particles may not pass typical track quality 
requirements and a dedicated analysis is required. 

While fractionally charged particles are common in some theoretical scenarios (e.g, super- 
strings 0131), a variety of searches for these objects in bulk matter, cosmic rays and accelerator 
based experiments have reported null results [4j. Strong constraints on models with fraction- 
ally charged states come from astrophysics and cosmology 0. These constraints, however, do 
not apply if the reheating temperature of the universe after the last stage of inflation is much 
lower than the mass of the fractionally charged particle such that there is no thermal relic from 
freeze-out [6j. In this Letter, we search for suchs particles, using as a benchmark the scenario 
considered in [5j, namely a model with new, fractionally charged, massive spin-1/2 particles 
that are neutral under SU(3)c and SU(2)i and therefore couple only to the photon and the Z. 
We denote such particles with fractional charge iqe, where e is the charge of the electron, as 
Lq, and assume they have a lifetime sufficiently long such that they do not decay within the 
detector volume. 

An interesting possibility is that these L q could also be charged under a new asymptotically- 
free gauge group SU(N) with a confinement scale A. This would make them a variant type of 
"quirk", which are quark-like, naturally fractionally charged particles [7j . For A < 100 eV, 
the effect of the confining string is invisible on collider detector scales and pair produced 
quirks would have the same kinematic properties and experimental signature as the bench- 
mark model considered in this Letter. However, the existence of this string would cause even- 
tual annihilation of any pairs formed in the early universe resulting in a negligible relic abun- 
dance [6 1, thereby evading the constraints cited by Ref. |5[ irrespective of cosmological history. 

We search for Lq, Lq particles in a sample of tracks with a muon-like signature. We identify 
fractionally charged particle candidates by their anomalously low ionization energy loss in the 
inner tracker. This study complements Compact Muon Solenoid (CMS) searches for heavy 
stable charged particles with anomalously high ionization energy loss fl8j|9l . 

The central feature of the CMS apparatus [ 10 1 is a superconducting solenoid of 6 m internal di- 
ameter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip 
tracker, a lead-tungstate crystal electromagnetic calorimeter, and a brass /scintillator hadron 
calorimeter. Muons are measured in gas-ionization detectors embedded in the steel return 
yoke. Extensive forward calorimetry complements the coverage provided by the barrel and 
endcap detectors. Of particular importance to this search is the inner tracker |TT][, which con- 
sists of 1440 silicon pixel and 15 148 silicon strip detector modules. The inner tracker records 
16 measurements on average per track. 

The analysis is performed on the pp collision data sample recorded by the CMS detector at 
y/s = 7TeV in 2011, corresponding to an integrated luminosity of 5.0 fb -1 . Pair production 
of LqLq at the LHC proceeds via a modified Drell-Yan mechanism with weak isospin t^i = 0, 
which has Lq-Z axial coupling g A = and vector coupling gy = —2q sin 2 Qyj. We have per- 
formed Monte Carlo simulations of this signal using PYTHIA v6.422 |T2|, with q = 1/3,2/3, 
and 1, and with masses of 100, 200, 300, 400, 500, and 600 GeV. The detector response is mod- 
eled with a simulation based on Geant4 [13J. 

The data are selected with a single-muon trigger that requires a track reconstructed in both the 
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inner tracker and muon detectors with transverse momentum pj > 40 GeV and pseudorapidity 
\t] | < 2.1, where tj = — ln[tan(0/2)] and 8 is the polar angle. The radius of curvature of a 
fractionally charged particle in a magnetic field is larger than that of a particle of unit charge 
with the same momentum, so the reconstructed pj is larger than the true pj by the inverse of 
the particle's charge. 

The trigger requirement that unit charge particles have pj > 40 GeV corresponds to a require- 
ment of pj > 27 GeV for L 2 /3 and pj > 14 GeV for I4/3. The trigger efficiency for L 2 /3 is in 
the range 67-74%, which is very similar to the efficiency for unit charge particles simulated 
with the same mass. By contrast, the I4/3 trigger efficiency is between 8% and 18%. The lower 
trigger efficiency for I4/3 results from the fact that many of the energy deposits in the tracker 
and muon detectors are below the threshold required to record a measurement. The trigger 
efficiency for I4/3 is a strong function of the particle's velocity. The reconstruction efficiency 
for slower-moving I4/3 particles is larger, since the energy deposits are more likely to be above 
threshold. As a result, the reconstructed velocity distribution is very different from the gen- 
erated distribution. For L 2 /3 the larger overall efficiency means that the velocity distribution 
is less affected by the reconstruction, but slower moving L 2 /3 particles fail the signal region 
requirement since their recorded dE/ dx measurements are too large. 

In a pre-selection step, candidates for fractionally charged particles are defined as tracks recon- 
structed in the inner tracker and matched to a track in the muon detectors. The pre-selection 
criteria, which are described below, are chosen to obtain well-reconstructed tracks and to sup- 
press background from cosmic ray muons. We consider candidate muon tracks with large 
reconstructed transverse momenta, pj > 45 GeV, as measured in the inner tracker, in the range 
\t] I < 1.5. A loose requirement on the track fit quality, ^ 2 /dof < 10, rejects very poorly re- 
constructed tracks. We also require at least six dE/ dx measurements from the tracker, where 
a dE/ dx measurement is the signal amplitude recorded in an inner detector module divided 
by the particle's path length through the module. To ensure that the track is isolated, the sum 
of the pj of all other tracks within a cone of AR = \J (A(p) 2 + (A//) 2 < 0.3, where <p is the az- 
imuthal angle, around the candidate must be less than 0.1 times the pj of the candidate. The 
sum of the electromagnetic and hadronic calorimeter energy recorded within this cone, includ- 
ing that deposited by the candidate, must be less than 5 GeV. Muons from cosmic rays are 
found in only a small fraction of the triggered events, but because they typically arrive out of 
coincidence with the bunch crossing, the charge is sampled away from the signal's maximum. 
The resulting low dE/dx measurement mimics those of a fractionally charged particle. Sev- 
eral requirements help to suppress the cosmic ray background. The primary vertex with the 
smallest distance to the point of closest approach of the candidate track is required to be a well 
reconstructed vertex [14J. The track must also have at least two dE/ dx measurements in the 
pixel tracker. The candidate track is required to satisfy \d z \ < 0.5 cm and \d xy | < 0.1 cm, where 
d z and d xy are the longitudinal and transverse impact parameters with respect to the primary 
vertex. The estimated time for the candidate particle to be at the interaction point, as measured 
by the muon detectors, under the assumptions that the particle has velocity /3 = 1 and is mov- 
ing from the interaction point outward, must not be earlier than the nominal bunch crossing 
that triggered the event. We require a max < 2.8 rad, where a max is the maximum 3-dimensional 
opening angle between the candidate track and any high-momentum (pj > 35 GeV) track re- 
constructed in either the inner tracker or in the muon system alone. 

The combined efficiency of the trigger and the pre-selection for signal events generated with 
mass 100 GeV is 45% for L2/3 and 4.4% for I4/3. Most of the efficiency loss is at the trigger stage, 
where a signal track may fail to be reconstructed. The pre-selection requirements reduce the 
signal efficiency by roughly a factor of two. 
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After the pre-selection, events with two Lq candidates are rejected from the search sample if the 
invariant mass of the two candidates nin is in the range 80 < mn < 100 GeV. This ensures that 
the search sample is independent of a control sample, defined below, that is used to estimate 
the collisions background. Although in the considered signal model, LqLq are produced in 
pairs, events containing a single candidate track are retained after the pre-selection. Events 
containing more than two candidates, which make up less than 0.1% of the search sample, are 
excluded. 

We isolate the signal using a technique that imposes few assumptions on any particular model 
but nonetheless has the power to suppress the large standard model backgrounds from pp col- 
lisions. A fractionally charged particle is most clearly distinguished from a standard model 
particle by its lower rate of energy loss in the detector since dE/dx <x q 2 [15 j. Figure [j] shows 
the distribution of dE/ dx measurements associated with tracks passing the pre-selection, for 
data and for the 100 GeV L2/3 and I4/3 simulated signal samples. The measured values from 
data lie predominantly in the region above 2 MeV/cm. We therefore define low-ionizing mea- 
surements to be those with dE/dx < 2MeV/cm. By requiring a large number of such low 
dE/ dx measurements, SM backgrounds can be suppressed while most of the signal events that 
pass the pre-selection are retained. 

A signal region is determined by maximizing the expected limit on the production cross section 
of fractionally charged particles while varying the minimum number of low dE/ dx measure- 
ments. The signal region optimization is performed for simulated samples with a Lq mass of 
100 GeV and 400 GeV, and for both samples the optimum signal region is defined with the re- 
quirement that events contain a track with at least six low dE/ dx measurements. Of the signal 
events passing the pre-selection, 75% of L2/3 events are in the signal region, and 93% of I4/3 
events are in the signal region. 

We use control samples of data to estimate the background contribution from cosmic ray muons 
and from particles produced in pp collisions. To estimate the cosmic ray background, we use 
a sample of muons obtained with the nominal pre-selection except for two inverted require- 
ments, 0.1 < \d X y\ < 1.1 cm and 0.5 < \d z \ < 50cm. The weight factor for this selection is the 
product of two weights, each determined from a cosmic ray enriched sample as the ratio of the 
yield in the nominal selection region to the yield after inverting a single requirement. The yield 
in the signal region is scaled to the expected yield for the nominal pre-selection. The cosmic 
ray background estimate is 0.007 events. 

A Z-peak control sample is used to estimate the pp collision background. This sample is se- 
lected by relaxing the transverse momentum requirement to pj > 35 GeV, requiring 80 < 
nin < 100 GeV, and applying all other pre-selection requirements. Figure [I] shows the dis- 
tributions of dE/ dx measurements associated with selected tracks for both the search sample 
and the Z-peak control sample data. The two distributions agree within the statistical uncer- 
tainties over the full dE/ dx range. The ratio of the number of tracks passing the pre-selection 
in the search sample to the number of tracks passing the pre-selection in the Z-peak control 
sample is 10.5. The control sample is scaled by this ratio to obtain the estimate of the collision 
background. 

To obtain a more precise estimate of the expected background in the signal region, we extrapo- 
late from the background-dominated region of events containing a pre-selected track with zero 
to five low dE/ dx measurements. For a muon from a Z decay, the dE/ dx measurements asso- 
ciated with the track are expected to be uncorrelated, and the number of measurements below 
a given dE/ dx value can be described by a generalized binomial distribution. The fit of this 
function to the background-dominated region of the Z-peak control sample is shown in Fig. [2] 
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Figure 1: Distribution of dE/ dx measurements associated with tracks passing the pre-selection 
in the search sample and the Z-peak control sample data. Simulated LqLq signal samples for a 
mass of 100 GeV are shown for q = 2/3, 1/3. The distributions are normalized to the area of 
the search sample. The magenta vertical line at dE/ dx = 2 MeV/cm indicates the upper limit 
of the range of measurements considered to be low-ionizing. 

Extrapolation of this function into the signal region yields a pp background estimate of 0.005 
events. 

The systematic uncertainties that significantly impact the results are the uncertainties in the 
integrated luminosity the background estimate, and the signal efficiency. The uncertainty in 
the integrated luminosity is 2.2% [Tl6| . 

The systematic uncertainty on the estimate of the background from cosmic rays is assessed by 
examining a sideband defined by replacing the inverted d z selection with a different inverted 
selection requirement, a max < 2.8 rad, and taking the difference with respect to the nominal 
determination as the systematic error, determined to be ±0.006 events. For the uncertainty in 
the pp background estimate, three sources are assessed: the statistical uncertainties in the fit- 
ted parameters of the binomial function, the choice of the function used to fit the background 
estimate in the background-dominated region, and the uncertainty in the model provided by 
the control sample to describe the background in the signal region. The three sources of un- 
certainty in the pp background estimate are summed in quadrature for a total uncertainty of 
±0.004 events. 

The systematic uncertainties in the signal efficiency are summarized in Table [l] The simulation 
of the trigger efficiency for a fractionally charged particle is sensitive to the accurate modeling 
of the muon detectors' electronics and gas gain as well as the threshold for recording a hit, 
because it has less ionization energy loss and thus the peak of its Landau distribution is closer 
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Figure 2: Number of events with at least one track with the given number of low dE/ dx mea- 
surements, for search data and the Z-peak control sample background estimate. The binomial 
function used to represent the background estimate is shown, with the band representing its 
uncertainty. 

to the threshold than that of a muon. We assess the impact of this uncertainty on the signal 
region efficiency by varying the gain in the muon system by a conservative estimate of its 
uncertainty. The impact of such a variation L 2 /3 is small, since the charge distributions are 
typically far above the threshold. However, for I4/3, the probability to record a hit degrades 
significantly as the gain decreases. The impact of this variation for L 2 /3 is 1%, and for I4/3 
is 8.5%. The systematic uncertainty in the offline global muon identification requirement is 
negligible by comparison. 

The scale of measured dE/ dx values impacts the track reconstruction efficiency, as well as the 
fraction of pre-selected tracks that pass the signal region selection. A discrepancy is observed 
between the dE/ dx distributions for simulated events and for events from the Z-peak control 
sample. The magnitude of the effect has been quantified by determining the change in scale of 
dE/ dx required to give the same fraction of events in the region below dE/ dx = 2MeV/cm 
in both samples. This scale change is found to be +5%, and to test its impact on the track re- 
construction efficiency, we vary the amplitude of the simulated dE/ dx measurements by ±5% 
before the track reconstruction is performed. A variation of +5% in the dE/ dx measurements 
changes the signal efficiency by +16% for I4/3 particles and —7.5% for L2/3. These variations in 
the signal efficiency are taken as the systematic uncertainties. The uncertainty resulting from 
the track momentum scale uncertainty is established to be less than 1%. The impact of the 
uncertainty in the muon timing measurements is assessed to be 2%. 

The numbers of expected background and observed events are summarized in Table |2j We 
observe zero events in the data search sample in the signal region, which is consistent with 
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Table 1: Systematic uncertainties (in %) in the signal efficiency. 



Source 


^1/3 


^2/3 


Muon trigger 


8.5 


1 


Track reco. (dE/ dx scale) 


16 


7.5 


Track momentum scale 


<1 


<1 


Muon timing 


2 


2 


Total 


18 


8 



Table 2: Numbers of events in the signal region, containing a track with at least six dE/dx 
measurements with dE/ dx < 2MeV/cm, estimated from the background and observed in the 
search data. 



Cosmic rays 


0.007 


± 


0.006 


pp collisions 


0.005 


± 


0.004 


Total background 


0.012 


± 


0.007 


Observed events 










the background estimate of 0.012 ± 0.007 events. The L q Lq signal efficiency and average ft for 
different signal masses and charges are listed in Table [3] Ninety-five percent confidence level 
(CL) upper limits on the L q L q production cross section are calculated using the CLs criterion 
[Tl7llT8| . Expected and observed 95% CL limits are shown in Fig. |3j We exclude the production 
of L 2 /3 with a mass below 390 GeV and the production of I4/3 with a mass below 250 GeV at 
95% CL. 

Table 3: The signal efficiency and average velocity (jS) for different mass points and charge 
hypotheses. 



mass 


L2/3 




^1/3 




(GeV) 


signal eff. 


(P) 


signal eff. 


(P) 


100 


0.341 ± 0.026 


0.84 


0.041 ± 0.007 


0.52 


200 


0.357 ± 0.027 


0.83 


0.060 ± 0.011 


0.51 


300 


0.337 ± 0.026 


0.82 


0.074 ± 0.013 


0.51 


400 


0.314 ± 0.024 


0.80 


0.091 ± 0.016 


0.51 


500 


0.265 ± 0.020 


0.79 


0.104 ± 0.019 


0.51 


600 


0.251 ± 0.019 


0.78 


0.109 ± 0.019 


0.51 



In summary using the signature of at least six low dE/ dx measurements, a search has been 
performed for heavy, long-lived, lepton-like fractionally charged particles. The search is based 
on a pp collision sample recorded by the CMS detector at y/s = 7 TeV, corresponding to an inte- 
grated luminosity of 5.0 fb _1 . The number of observed events in the signal region is consistent 
with the background estimate. Upper limits on the production cross section of pair produced, 
spin-1/2 particles that are neutral under SU(3)c and SU(2)i exclude at 95% confidence nil 
below 390 GeV for particles with q = 2/3 and mi below 250 GeV for particles with q = 1/3. 
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Figure 3: Expected and observed limits on the cross section of L^L^ for q = 2/3 and 1/3. The 
theoretical prediction of modified Drell-Yan production with t$i = is shown. The lines for the 
expected and observed limits are overlapping. 
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